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ABSTRACT 

Aims. We present the first magnetic Doppler images of a rapidly oscillating Ap (roAp) star. 

Methods. We deduce information about magnetic field geometry and abundance distributions of a number of chemical elements on 
the surface of the hitherto best studied roAp star, HD 24712, using the magnetic Doppler imaging (MDI) code, INVERSIO, which 
allows us to reconstruct simultaneously and consistently the magnetic field geometry and elemental abundance distributions on a 
stellar surface. For this purpose we analyse time series spectra obtained in Stokes I and V parameters with the SOFIN polarimeter 
at the Nordic Optical Telescope and recover surface abundance structures of sixteen different chemical elements, respectively ions, 
including Mg, Ca, Sc, Ti, Cr, Fe, Co, Ni, Y, La, Ce, Pr, Nd, Gd, Tb, and Dy. For the rare earth elements (REE) Pr and Nd separate 
maps were obtained using lines of the first and the second ionization stage. 

Results. We find and confirm a clear dipolar structure of the surface magnetic field and an unexpected correlation of elemental 
abundance with respect to this field: one group of elements accumulates solely where the positive magnetic pole is visible, whereas 
the other group avoids this region and is enhanced where the magnetic equatorial region dominates the visible stellar surface. We also 
observe relative shifts of abundance enhancement- or depletion regions between the various elements exhibiting otherwise similar 
behaviour. 
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1. Introduction 

Many physical processes important for stellar evolution are cru- 
cially influenced by stellar magnetic fields. Theoretical and ob- 
servational frameworks tell us for instance that they influence 
microscopic diffusion of chemical elements within stellar atmo- 
' spheres, and they may also redistribute angular momentum or 
give rise to enhanced hydrodynamical instabilities and thus mod- 
ify mixing properties in stellar interiors. 

Ap stars, representing about 5-10% of the upper main se- 
quence stars, exhibit magnetic fields that appear to be highly 
ordered, very stable and often very strong. Many Ap stars also 
show strong spectral line profile variations synchronized to stel- 
lar rotation. This is attributed to oblique magnetic and rotation 
axes and to the presence of a non-uniform distribution of chemi- 
cal elements on their surface within the framework of the oblique 
rotator model (introduced by Stibbs, 119501 ). With a few excep- 
tions such inhomogeneities exist only in the atmospheres of A 
stars with magnetic fields, demonstrating that these fields play a 
crucial role in their formation and evolution. 

The spectra of Ap stars also exhibit a remarkable variety of 
sometimes unidentified spectral line features. Ryabchikova et al. 
( |2004| l for example find overabundances of some Fe peak el- 
ements, which may reach a few dex for Cr and even higher for 
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the REE, especially for their second ionization stage, while other 
chemical elements are found to be underabundant compared to 
the solar value. An important subgroup of the Ap stars, the 
rapidly oscillating Ap stars, in addition exhibits high-overtone, 
low-degree, non-radial p-mode pulsations with periods of 6-2 1 
minutes (Kurtz & Martinez 120001) . Their observed pulsational 
amplitudes are modulated with the magnetic field variations, ex- 
hibiting maximum amplitude at the phases of highest magnetic 
field strength (Kurtz , 1982J . This indicates a close connection be- 
tween the magnetic field and the excitation mechanism of roAp 
pulsations. 

The chemical peculiarities mentioned are attributed to the se- 
lective diffusion of ions under the competitive action of radiative 
acceleration and gravitational settling (Michaud 1970 ) possibly 
in combination with a weak, magnetically driven wind (see, e.g., 
Babel l 19921 1. 

We focus our research on roAp stars since these objects ex- 
hibit a number of new, largely unexplored phenomena related 
to their pulsational variability while sharing all the properties 
of other Ap stars. The roAp stars are the only main-sequence 
stars in which high-overtone p-mode pulsations are easily de- 
tected and pulsation modes can be identified. Furthermore, roAp 
stars represent the only group of pulsating stars for which mag- 
netic and pulsation characteristics can well be constricted ob- 
servationally. Hence, peculiar pulsating stars are key objects in 
the investigation of the atmospheric and internal structure of the 
middle main sequence stars. 

Very little is known as yet about the origin and structure of 
magnetic fields and their connection and interaction with surface 
abundance patches and pulsation. With the development of high- 
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resolution spectropolarimeters it has become possible to extract 
the full amount of information about magnetic field and abun- 
dance distributions using Stokes parameter observations and ap- 
plying magnetic Doppler imaging. 

We have developed and tested a code for MDI 
(Piskunov&KochukhovEOOll Kochukhov & Piskunov |2002ll 
that allows us to reconstruct simultaneously and consistently 
the magnetic field vector and abundance distribution of various 
elements on a stellar surface by taking into account all relevant 
physics of polarized line formation. 

In deriving the surface abundances of numerous light-, iron- 
peak- and rare earth elements by applying MDI, we get more in- 
sight into the relation between the magnetic fields, vertical and 
horizontal abundance characteristics and pulsation, and thus into 
the atmospheric structure of Ap and roAp stars. At present sur- 
face abundance structures have been derived for only one other 
roAp star, HD 83368 (Kochukhov et al. 2004b). 

Our paper is composed in the following way: Sect.|2] is de- 
voted to the description of the target star, HD 24712, on which 
we comment in Sect.[3]regarding our spectropolarimetric obser- 
vations and the data reduction. Details about MDI, the magnetic 
field geometry of HD 24712, and the surface abundance struc- 
tures of individual elements can be found in Sect.|4]while Sect.|5] 
is devoted to the discussion of our results. 

2. HD24712 

HD 24712 (HD1217, DO Eri) is o ne of t he best-studied roAp 
stars with light (Wolff & Morrison 119731 . spectrum and mag- 
netic variations that was discovered to pulsate by Kurtz ( 1982). 
Matthews et al. (.1988) found radial velocity variations with an 
amplitude of 0.4+0.05 km s 'and the main photometric period 
of 6.14min. 

During an observing campaign with the Whole Earth 
Telescope (WET), Kurtz et al. ( 120021 ) detected a 'missing' mode. 
The detailed investigations by Cunha et al. (I2003l and l2006 ) sug- 
gest in this context, that the p-mode pulsations are strongly af- 
fected by the global stellar magnetic field. 

A study of the surface abundance inhomogeneity in 
HD 24712 based on the observed line profile variations was 
performed by Ryabchikova et al. ( I2000I I. Surface element dis- 
tribution was schematically approximated as a set of circular 
spots whose position, radii and the abundances inside and out- 
side the spots have been obtained with an iterative procedure. 
Ryabchikova et al. (I2000I I found that Pr ii, iii, Nd ii, iii, and Co 
are concentrated in one large spot near the magnetic pole, while 
Fe is concentrated rather in 4 spots around the magnetic equator 
They also found that centers of Nd, Pr and Co spots do not co- 
incide with the pole of the magnetic field geometry determined 
by Bagnulo et al. (1995), but are shifted by » 10° towards the 
magnetic equator. According to their analysis Pr-Nd and Co are 
overabundant by a; -Hi dex in the spot relative to the outside sur- 
face area and one needs x2 dex abundance differences between 
Pr and Nd in the first and the second ionization stages. The max- 
imum surface abundance gradient for Fe was 0.6 dex. 

We performed a detailed spectroscopic pulsation study of ra- 
dial velocity variations in high resolution and high time-resolved 
data (obtained with SOFIN at NOT and with UVES at the ESO 
VLT). Simultaneously we performed high precision photomet- 
ric observations from space with the MOST (Microvariability 
and Oscillation of Stars) satellite (Walker et al. |2003i ) During 
this process, we discovered a detailed pulsation picture of the 
star (Ryabchikova et al. 120071) . We found maximum radial ve- 
locity amplitudes for rare earth lines and in addition phase shifts 



Table 1. Basic physical parameters of HD 24712. 



T^ff 7250+ 150 a: 

logg 4.2 + 0.1 

M 1.63 Mo 

L/Lq 0.91 ± 0.04 

ProtCd) 12.45877 ± (16) 

^'puis ==6.15 min 

VeSin; 5.6 km s"' 



progressively increasing from one element to another and within 
different groups of lines of the same element. We propose that 
this pulsational behaviour is caused by the fact that these ele- 
ments are concentrated in a thin layer in the stellar atmosphere 
where pulsation reaches its maximum amplitude. 

These observations and a possible correlation to theoretical 
investigations made HD24712 an ideal candidate for spectro- 
scopic and spectropolarimetric observations with ground based 
equipment. 

The basic physical parameters of HD24712 are listed in 
Table[T] T^g and log g were applied as derived by Ryabchikova 
et al. (1 19971 ). The atomic data used in our study was extracted 
from the Vienna Atomic Line Database (VALD, Piskunov et al. 
[T9951 Ryabchikova et al. [19991 Kupka et al. 17999] ). 

3. Observations and data reduction 

Spectropolarimetric observations of HD24712 were carried out 
in 2003, between Oct. 29 and Nov. 18, using the high-resolution 
echelle spectrograph SOFIN attached to the Cassegrain focus of 
the 2.56 m Nordic Optical Telescope (NOT), La Palma, Spain. 
The spectrograph is equipped with three different cameras, and 
to obtain observations in the spectropolarimetric mode, the sec- 
ond camera was used. For MDI we used 49 different spectral 
lines from seven spectral orders, each covering about 40 A to 
50 A between 4000 and 6000 A. 

The polarized spectra were obtained with a Stokesmeter, 
consisting of a calcite plate used as a beam splitter and an achro- 
matic rotating quarter-wave plate. The position of the quarter- 
wave plate is controlled by a stepping motor. Usually a sequence 
of four subexposures is processed in order to obtain accurate 
circular polarization measurements, where each of the beams is 
exposed twice, with the quarter-wave plate rotated by 90° after 
the first and before the last exposure. This way it is possible to 
reduce instrumental effects to a minimum, because in the images 
taken with the quarter-wave plate rotated by 90°, possible instru- 
mental signatures change sign and cancel each other out during 
the averaging of the two exposures. 

For data reduction, the 4A software package written in C 
(Ilyin 2000) was used. All standard procedures like bias sub- 
traction, flat field correction, subtraction of the scattered light, 
weighted extraction of the orders, and bad pixel (cosmic ray) 
corrections were included. ThAr exposures obtained before and 
after each observing night were used to perform wavelength cal- 
ibration and to test for possible spurious instrumental polariza- 
tion, caused for instance by bending of the spectrograph, which 
is directly mounted on the telescope, different positions of the 
star on the slit or temporal variations of the seeing. S/N ratios 
for the observed spectra typically reached 500-600 per pixel 
at A- 5000 A, except for one spectrum which has S/N 300. 
Reduced spectra have a resolving power of ^ 40 000. 
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Table 2. Journal of spectropolarimetric observations of 
HD 24712. 



HJD Rotational phases 5 IN 

(245 0000+) 



2941.6516 


0.4396 


290 


2943.6341 


0.6001 


550 


2945.5758 


0.7562 


520 


2946.5973 


0.8380 


510 


2947.6067 


0.9171 


570 


2948.6514 


0.0026 


660 


2952.6267 


0.3207 


560 


2953.6572 


0.4040 


520 


2954.6316 


0.4824 


550 


2955.6350 


0.5620 


570 


2956.5954 


0.6390 


620 


2957.5958 


0.7204 


560 


2961.6160 


0.0414 


500 



In total we obtained observations at 13 rotational phases. 
Each pair of the Stokes / and V spectra was calculated from a 
sequence of polarimetric exposures obtained over 35 min. This 
is sufficiently long to average out variations due to the 6.15 min 
pulsation period of HD 24712. 

Rotation phases of HD24712 (see Table |2]l were calculated 
according to the ephemeris and rotation period derived by 
Ryabchikova et al. (|2005T l: 

HJD«BJ„„) = 2453235.18(40)+ 12.45877(16)d. 

4. Magnetic Doppler imaging (IVIDI) 

Doppler imaging (DI) has become a very successful tool to gain 
information on inhomogeneities on the surface of a rotating star 
caused by changing abundance and/or magnetic field structure. 
The first method for reconstruction of stellar magnetic field ge- 
ometries from Stokes V and I spectropolarimetry was developed 
by Piskunov (1985) and applied to Ap stars by Glagolevskii et al. 
(1985). A similar technique has been very successfully applied 
to mapping magnetic field topologies in late-type active stars by 
Donati & ColHer Cameron (1997), Hussain et al. (2001), and 
Donati et al. (2006). 

By inverting time series of spectropolarimetric observations 
of rotating stars based on regularized image-reconstruction pro- 
cedures implemented in IN VERS 10 developed by Piskunov 
& Kochukhov (2002), it became possible to invert simulta- 
neously the information of the rotationally modulated Stokes 
profiles into elemental abundance as well as the three vector 
components of the stellar magnetic field. In their papers I and 
II, Piskunov & Kochukhov (2002) and Kochukhov & Piskunov 
( I2OO2I 1 describe in detail the concept of MDI, the basic numerical 
techniques to efficiently solve the magnetic radiative transfer, the 
surface integration of ffie resulting Stokes profiles and the appli- 
cation of the regularization concept. Applications of INVERS 10 
are presented in Kochukhov et al. (2002, 2004a) and Folsom et 
al. (l2008] l. 

In case of HD 24712 we only have access to two Stokes pa- 
rameters (I and V). For ffiis reason we had to resort to multi- 
polar regularization to obtain a solution free from undesirable 
cross-talks between different field components. This is a rea- 
sonable prior model of the Ap-star magnetic field geometry. 
Numerous previous studies (e.g. Landstreet & Mathys 120001 



Bagnulo et al. |2002| l showed that rotational modulation of the 
integral quantities derived from the circular polarization spec- 
tra of Ap stars are well-explained by dipolar or quadrupolar- 
like fields with some additional influence of abundance inhomo- 
geneities. Multipolar regularization (see Piskunov & Kochukhov 
I2OO2I 1 has been successfully used for the mapping of the Ap 
stars g^CV n (Kochukhov et al. 2002) and HD 72106 (Folsom 
et al. I2OO8I 1. This method introduces external constraints simi- 
lar to the multipolar field parameterization proposed by Bagnulo 
et al. d 19991 1 and leads INVERS 10 to solutions close to a 
general second-order multipolar expansion. Numerical exper- 
iments (Kochukhov & Piskunov 2002 ) demonstrated that this 
technique provides reliable imaging of stars with quasi-dipolar 
or quadrupolar magnetic fields. However, our previous experi- 
ence with the MDI shows that using multipolar regularization 
may prevent us from reconstructing the fine structure of mag- 
netic field by forcing it closer to a low-order multipolar geom- 
etry. The alternative, Tikhonov regularization, which does not 
introduce assumptions about the global geometry of ffie stellar 
surface structures, was only applied for reconstruction of chemi- 
cal abundance. The values of two regularization parameters were 
adjusted to ensure that at the convergence the total contribution 
of the regularization terms was roughly 20-80% of the ffie sum 
of weighted squared differences between observations and syn- 
thetic spectra. We have verified that a variation of regularization 
parameters by factors of 2-5 leaves little imprint on the final so- 
lution. 

The additional information content of the Stokes profiles 
compared to the data used for conventional Doppler Imaging 
makes it possible to also map stars with rotational velocities 
as low as that of HD 24712 with Vg sin /=5.6 km s"' (derived by 
Ryabchikova et al. I1997I I. Nevertheless we would like to state 
here that due to this very low Vg sin /, the uncertainties of the 
latitudes of the abundance enhancement and depletion regions 
are larger than for stars with Vg sin / higher than ^10 km s ' . On 
the other hand, the longitude position of the surface features is 
mainly constrained by the phase coverage, which is fairly good 
for our observations. Information contained in the considerable 
rotational modulation of the line intensities due to changing vis- 
ibility of abundance spots on the surface of HD 247 12 allows us 
to infer the longitudinal position of large-scale abundance struc- 
ffires wiffi sufficient precision. Numerical experiments demon- 
strate that longitudinal shifts of major abundance features by 
~ 10° (few 0.01 in phase units) can be readily detected with our 
observational data. 

4. 1 . Magnetic field geometry 

Determining ffie geometry of ffie magnetic field on the surface of 
HD 24712 was performed choosing seven different Fei and five 
different Nd m lines, which are listed in Table[3] These lines were 
chosen because their Lande factors and Zeeman splitting pat- 
terns yield strong polarization signatures, they are not affected 
by severe blending, and they show clear variations in the strength 
of ffie Stokes I profiles due to abundance inhomogeneities and, 
to a lesser extent, due to variable Zeeman effect. Furthermore, 
the lines of Fe and Nd exhibit an opposite variation with rota- 
tion phase. Thus, using these two elements simultaneously in a 
self-consistent abundance and magnetic inversion yields a more 
robust reconstruction of the stellar magnetic field topology. 

Important parameters needed as input for the MDI proce- 
dure are the inclination and the azimuth angle of the stellar rota- 
tional axis, / and 0. It should be mentioned that when mapping 
with Stokes I and V components, the azimuth angle is irrelevant 
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Fig. 1. Mapping of the distribution of magnetic field strength (a) and field orientation (b) on the surface of the roAp star HD 24712, 
as derived with modelling the Stokes I and V parameter stellar observations. The results reveal a dominant dipolar magnetic field 
geometry on the surface of this star with the field strength varying between 2.1 kG and 4.2 kG. 



ip=(i.m ^=0.20 p=0.40 p=0.f/Q ^1=0.80 




Fig. 2. The abundance distributions of Fe and Nd on the surface of HD 24712. The upper panel shows the surface map of the 
Fe abundance and the lower panel presents the distribution of Nd derived from Ndiii lines. These maps were derived using the 
Stokes I and V parameter spectra as presented in Fig. [3] In this and other similar plots we use the s(Xi) to denote abundance, in the 
log(A^x/Mot) units of the element "X", derived from lines belonging to the ionization stage "i". 



and we cannot distinguish between / and 180° - i. However, we 
adopted / > 90° to be consistent with the analysis of Bagnulo et 
al. (1995) , who determined both angles based on the traditional 
circular- and on broadband linear polarization observations and 
found / - 137° and = 4°. Adopting these values, a sur- 
face magnetic map very close to a dipolar geometry yielded the 
best fit to the observed line profiles. The magnetic field strength 
varies between 2.1 kG and 4.2kG. In Figs. [1] and [3] the observed 
line profiles of the chosen Fei and Ndiii lines (dots), the cor- 
responding fit reached during the inversions (lines) and the re- 
sulting magnetic field geometry are presented. The near dipolar 
structure agrees very well with the dipolar model proposed by 
Bagnulo et al. ( |1995l l, only differing in the polar field strength, 
where they obtain Bp = 3.9 kG, compared to 4.2 kG derived from 
magnetic Doppler Imaging. 



4.2. Fe and Nd abundance distribution 

The surface abundances of Fe and Nd mapped simultaneously 
with the stellar magnetic field geometry are shown in Fig.|2] 
As mentioned above, these distributions were calculated using 
seven Fei and five Nd iii lines listed in Table[3] One can easily 



see that abundances of both elements are quite globally struc- 
tured. Fe varies between -5.3 and -4.7 dex in the log(A^/A^tot), in 
its maximum, still slightly underabundant compared to the solar 
value (-4.59 dex). Nd abundance being derived from Ndm lines 
is very high and it varies between -8.0 and -6.9 dex (compared 
to -10.59 dex solar value). As we see further in Sect. 14. 3.31 Nd 
abundance distribution from Nd ii lines is similar in shape but by 
1 .5 dex smaller in magnitude following the known REE anomaly 
observed in all roAp stars (Ryabchikova et al. 120011 12004I I. Fe 
and Nd seem to be perfectly anticorrelated: Fe is accumulated 
where Nd is depleted, and minimum Fe abundance can be found 
where Nd is at its maximum. 



Comparing the patterns of these two elements to the derived 
magnetic field geometry presented in Fig.[T] we find that the Fe 
abundance enhancement region is associated with the area where 
the magnetic equator dominates the visible stellar hemisphere, 
around ^ 0.5, whereas the Nd map exhibits its area of max- 
imum abundance around the phase where the positive magnetic 
pole is oriented towards the observer Due to the orientation of 
the star we are not able to directly observe the negative magnetic 
pole. 
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It is surprising to see that, contrary to our expectations, both 
elements are not accumulated or depleted on both magnetic 
poles or around the magnetic equator. Furthermore, the center 
of maximum or minimum abundance, respectively, seems to be 
shifted in latitude and/or longitude with respect to the magnetic 
polar regions, not only for Nd, but also for several other elements 
presented in Sect. |4.3l 

4.3. Surface abundances of individual elements 

After having derived the magnetic field geometry of HD 24712 
using the Fei and Ndiii lines, we were able to recover abun- 
dances of sixteen additional ions applying the IN VERS 10 code 
and the previously determined fixed magnetic map. 

4.3.1. Magnesium, calcium, and scandium 

The surface abundance distribution of magnesium was modelled 
using one Mgi line at 5528 A with IN VERS 10, and we found a 
remarkable variation between -6.0 and -4.5 dex. Mg is, as Sc, 
Ti, Fe, and Ni depleted (underabundant by up to 1 .5 dex com- 
pared to the solar value) at the phase where the positive mag- 
netic pole dominates the visible stellar surface and shows a sin- 
gle, huge overabundance region around the magnetic equatorial 
zone, indicating that the center of maximum abundance does not 
coincide with the (suspected) negative magnetic pole. 

For calcium as well as for chromium (see Section l4.3.2l ) we 
found only a marginal abundance variation not exceeding 0.2 
dex and consequently do not consider our mapping results as 
significant. 

Scandium exhibits a variation of 1.2 dex over the stellar sur- 
face, between -8.9 and -10.1 dex, also accumulated where the 
magnetic equatorial region dominates the visible stellar surface, 
whereby the central region of the depleted area occurs at similar 
longitude as that of Mg and Fe and the central region of maxi- 
mum abundance precedes that of Mg and Fe by about 60°. The 
region of enhanced Sc abundance seems to occupy a larger frac- 
tion on the surface of HD 24712 than numerous other elements. 
At its maximum abundance phase, the Sc value corresponds to 
the solar one of 8.99 dex. 

4.3.2. Iron-peak elements 

Two lines of ionized titanium were used to recover the surface 
abundance structure of this element (Fig. |4|. As Mg and Sc, it 
shows a tendency to accumulate near the magnetic equatorial 
region, exhibiting near solar concentration, and the abundance 
decreases to -7.8 dex where the positive magnetic pole is visible. 

For chromium only one line around 5510 A could be used as 
input for INVERSIO and we obtained just a marginal (insignif- 
icant as for calcium, see Section |4. 3. 2b variation between -5.9 
and -6.1 dex. 

Cobalt, as Y and all the REE, is clearly accumulated on the 
surface region visible at the phase where the positive magnetic 
pole dominates, and less abundant around the magnetic equato- 
rial region, varying between -5.8 to -6.4 dex. Co is significantly 
overabundant with respect to the Sun (-7. 12 dex) over the whole 
stellar surface. The element was mapped taking into account six 
Coi lines shown in Fig. |4] and listed in Table[3] 

A very interesting example of enormous line profile varia- 
tion has been discovered for the nickel lines. Nil at 5035 A, and 
5 146 A. The surface map derived from the single unblended line 
at 5035 A reveals a huge, high contrast region of overabundance 



at phase 0.5 and a region of underabundance down to -7.6 dex at 
the phase where the positive magnetic pole dominates the visible 
stellar surface. We noticed that the Nil /15 146.48 A line totally 
disappears at this rotational phase. In order to check the abun- 
dance derived from the single 5035 A line and to be able to use 
the blend of Ni with Fe, Co, and Ce around 5 146 A, we used 
the abundance pattern obtained from the Nil 5035 line to map 
Co and Ce from the blend mentioned above and got a very good 
agreement of observed and calculated profiles. The fit to the line 
profiles reached from this inversion is presented in Fig. ID Ni dis- 
tribution, in longitude and in latitude, is very similar to that of 
Fe, and abundance varies from -7.6 to -6.3 dex, being signifi- 
cantly underabundant compared to the solar value of -5.81 dex. 

4.3.3. Yttrium and the rare-earth elements 

Fffr/wm varies from -9.0to -8.5 dex (overabundant up to 1.3dex 
compared to the solar value of -9.83) and exhibits a structure 
closer to that of the REE, with the center of abundance maxi- 
mum shifted to the northernmost part of the visible stellar hemi- 
sphere. Yttrium reaches its maximum and minimum, respec- 
tively, slightly after that of Nd. As presented in Figs. |6] |2l and 
[8] all the REE we mapped are enhanced around the phase where 
the positive magnetic pole is visible and depleted where the mag- 
netic equatorial region dominates the surface. For REE depletion 
means the minimum of the observed surface abundance which 
usually exceeds the solar photospheric abundance. We obtained 
abundance distributions of La, Ce, Pr, Nd, Gd, Tb, and Dy. For Pr 
and Nd, abundance distributions were calculated from the sep- 
arate analysis of the lines of the singly and doubly ionized ele- 
ment. 

Lanthanum is significantly overabundant compared to the so- 
lar value (-10.91) and varies from -9.8 to -8.8 dex. It covers 
stellar latitudes very similar to those of Co, Nd, and Gd, although 
due to the very low Vg sin / of the star, the latitudes of the abun- 
dance enhancement and depletion regions are affected by larger 
uncertainties than for stars with Vg sin / higher than ^ 10km s ' . 
The center of maximum abundance slightly lags behind that of 
Nd. 

Variation of Y and La lines is illustrated in Fig.|4] 
The cerium surface map was reconstructed using two Ce ii 
lines, also exhibiting the typical overabundance of REE in roAp 
stars of more than a dex, but varying only moderately over the 
surface between -9.3 and -8.9 dex. 

Two independent maps of praseodymium were derived using 
the lines of the singly and doubly ionized element. It is caused 
by the known REE anomaly observed in roAp stars, where av- 
erage abundance derived from these lines differs by 1.5 - 2.0 
dex (Ryabchikova et al. 2004), and it is not possible therefore to 
fit line profiles of both ionization stages within the same inver- 
sion. Prii /15136 A and two spectral lines of Priii around 49 10 A 
and 4929 A were used for inversions. Pr varies moderately be- 
tween -10.0 and -9.6 dex being derived from one Prii line, 
while higher abundances but with the same distribution shape 
are obtained from Priii lines. The Pr spot seems to be shifted in 
latitude relative to the Nd enhancement region. 

Abundances for neodymium were also derived separately 
from singly and doubly ionized lines due to the large abundance 
difference between the two ionization stages. The Nd surface 
map was derived from eight lines of Nd ii listed in Table[3] and 
exhibits the same distribution as that obtained from Nd iii lines 
(see Sect l4.2| ). Nd abundance varies from -9.3 to -8.5 dex. As 
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Fig. 3. Comparison of the observed (dots) and calculated (lines) Stokes I and V profiles. Spectra are shown with increasing phase 
value (listed in Table |2]i from top to bottom. The bars at the lower left and middle show the vertical (5%) and horizontal (0.5 A) 
scale. 



for the Pr case, the abundance difference derived from Nd ii and ter of maximum abundance shifted by ^ 30° to higher phases 

Nd m lines reaches 1 .5 dex. compared to Co, Nd, and Dy. 

Gadolinium is concentrated in a spot very similar to that of Dysprosium can also be found in a spot near the positive 

Nd, close to the positive magnetic pole region, varying moder- magnetic pole with the center of abundance slightly preced- 

ately between -8.5 and -8.8 dex, and is overabundant nearly up ing that of e.g. Gd. The element is varying between -8.3 and 

to 2.5 dex relative to the solar value of -10.92 dex. -8-8 dex and found to be overabundant by more than 2.0 dex 

The doubly ionized terbium line at 5505 A was used to map relative to the solar value of - 10.90 dex. 
the surface distribution of this element. Tb shows a strong over- 
abundance of up to -7.7 dex (compared to the solar value of g DisCUSSion 
-11. 76 dex) on the surface part appearing at the phases where 

the positive magnetic pole dominates the visible stellar hemi- Observing line profiles in several polarization states of stellar ra- 

sphere. It is modulated between -9. 1 and -7.7 dex, with the cen- diation and modelling their shapes with magnetic Doppler imag- 
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Fig. 4. Comparison of the observed (dots) and calculated (lines) Stokes I and V profiles used for mapping of Ti, Co, 
Spectra are plot with increasing phase value (listed in Table|2]| from top to bottom. The bars at the lower left and middle 
the vertical (5%) and horizontal (0.5 A) scale. Note the total disappearance and reappearance of the Nil line at 5146 A. 



Y and La. 
again show 



ing, we could derive the magnetic field geometry and abun- 
dance distributions of numerous elements inhomogeneously dis- 
tributed over the surface of the roAp star HD 24712. This is the 
first analysis of this kind for a magnetic pulsating star. Our inves- 
tigation shows a complex picture of the interplay between mag- 
netic field, pulsation and atmospheric inhomogeneities in this 
prototypical roAp star. However, the low Vg sin / of the star and 
the limited spectral resolution of our spectropolarimetric obser- 
vations did not allow us to reliably recover all but the largest 
spatial scales of the abundance and magnetic structures. For this 
reason we will discuss only the gross properties of the chemical 
spots and the magnetic field, without touching upon the higher- 
order moments of the surface structures. 

We obtain a dominant dipolar magnetic field structure 
(which might be more complex, though the general picture is 
certainly dipolar) and elemental abundance patterns that are cor- 
related to this geometry in an unexpected way: instead of abun- 
dance enhancement regions on both magnetic poles or around 
the magnetic equator, we observe huge enhancement or deple- 



tion regions around either the phase where the positive magnetic 
pole is visible or where the magnetic equatorial region domi- 
nates the visible stellar surface. 

In addition, even more exciting and novel is the fact that 
these enhancement or depletion regions for various chemical el- 
ements are shifted in longitude relative to each other. The center 
of Mg underabundance as well as that of Sc, Fe, and Ni can be 
found at phase 0.0, while that of Ti occurs around phase 0.082. 

The center of maximum abundance of Tb appears very close 
to phase if - 0.0, that of Co and Gd slightly after ip - 0.0, around 
0.02, Nd spots appear at (/? = 0.04, Pr and Dy at around ip = 0.06. 
Y and La are already shifted by 30° or » 0.082 in the direction 
of increasing phase value. Ce, whose center of maximum abun- 
dance lags behind all other elements, is shifted by as much as 
a; 0.12 or «i 45°. 

To investigate the reliability of these longitudinal offsets of 
the abundance structures we have performed a forward calcula- 
tion for Y, La, and Ce showing most of the conspicuous shifts. 
Using the final abundance maps reconstructed for these elements 
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we have synthesized their Hne profiles, rotating surface struc- 
tures in the longitudinal direction by small amounts and com- 
pared the resulting modification of Stokes / with the typical ob- 
servational uncertainty. This comparison is illustrated in Fig. |9] 
(online material). As is evident from this numerical experiment, 
the high S /N of our observations enables us to easily detect 
shifts of <: 0.05 for Y and Ce and even smaller shifts for La. Thus 
the relative longitudinal positions of spots of different elements 
are precise to within 0.02-0.05 of the rotation period. 

The diverse surface abundance behaviour observed in 
HD 24712 most likely gives us information on how the magnetic 
field influences the chemical diffusion of the various species 
within the stellar atmosphere. Accordingly we tried to correlate 
our results with simple atomic parameters (e.g., atomic weight 
or excitation potential), but did not find any clear evidence for 
such a connection. This gives rise to the speculation, as the ef- 
fect is clearly present, that there is a more complex correlation of 
the various atomic characteristics influencing chemical diffusion 
through stellar atmospheres in presence of a magnetic field. Thus 
these relative shifts in longitude of the mapped elements provide 
valuable constraints for diffusion models that try to explain how 
the magnetic field and other hydrodynamical processes regulate 
atmospheric diffusion. 

Moreover, we obtain shifts in latitude of abundance en- 
hancement regions exhibiting otherwise comparable accumula- 
tion tendencies. Sc, Ti, Fe, and Ni in this context appear in simi- 
lar latitudes, whereas Mg and Co seem to be shifted to lower lat- 



itudes by 30° and 45°, respectively, compared to the above men- 
tioned elements. Among the elements belonging to the group 
that is enhanced around the positive magnetic pole, Tb and Y 
seem to cover the highest latitudes, Ce and Pr are found to be 
shifted to lower latitudes by around 15°, and La, Nd, Gd, and Dy 
spots are shifted even more southwards by x 30°. 

This finding helps to explain the phase shift between pul- 
sation curves of spectral lines belonging to different REE el- 
ements, for instance, Pr and Nd, that we found in the recent 
study by Ryabchikova et al. ( I2007I I. According to the NLTE 
analysis of Nd ii, iii and Pr ii, iii line formation in the atmo- 
sphere of HD 24712 (Mashonkina et al. 2005, 2009), the lines 
of both elements are formed at nearly the same atmospheric 
layers. However, the pulsation analysis of Ryabchikova et al. 
(120071) shows that they are shifted by 0. 1 in the pulsation phase 
relative to each other. According to the theoretical calcula- 
tions by Sousa&Cunha (2008) and numerical simulations by 
Khomenko & Kochukhov (2009), the amplitude and phase dis- 
tributions with height in roAp atmospheres should depend on 
the local magnetic field strength and inclination. Slightly differ- 
ent location of Pr and Nd spots relative to the magnetic field 
pole may therefore be a reason for the pulsation amplitude and 
phase difference observed in Pr and Nd lines. This illustrates the 
potential and importance of MDI when interpreting pulsational 
observations of roAp stars. 

Vertically stratified abundances are very likely present in the 
atmosphere of HD 24712 (e.g., Ryabchikova et al. |2008l ). In par- 
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Fig. 6. Surface abundance distribution of Ni, Y, La, and Ce (from top to bottom) on the surface of HD 24712. The maps were 
obtained taking into account the Unes Usted in Table |3] 
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Fig. 7. Elemental abundance distribution on the surface of HD 247 12. Maps of Pr and Nd derived separately from the lines of the 
first and second ionization stages of both elements are presented. 
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Fig. 8. Surface abundance distributions of Gd, Tb, and Dy derived from the lines listed in Table [3] 



ticular, a major discrepancy between two ionization states of the 
same REE element is attributed to stratification (Mashonkina et 
al.i2009). To estimate the impact of stratification, we performed 
an analysis of the vertical gradient of Fe abundance (Liiftinger 
et al. 1200 4^ during the phase where the magnetic equatorial re- 
gion dominates the stellar hemisphere and also where the pos- 
itive magnetic pole is visible, choosing Fe lines that are most 
sensitive to stratification. Fe turned out to be significantly strat- 
ified, with abundance changing by * 2 dex between the upper 
and lower parts of the atmosphere of HD 24712. We also found 
a slight difference in the location of the abundance jump compar- 
ing results from the two extreme phases. In order to keep MDI 
results influenced as little as possible by the stratification effects, 
we chose lines for mapping that are hardly affected by a chang- 
ing vertical distribution. It is still possible however that part of 
the horizontal abundance structure we find is due to variation of 
the chemical stratification profile across the stellar surface. 

Currently, there are no theoretical diffusion models that can 
provide quantitative explanation of the abundance maps derived 
in our Doppler imaging study of HD 24712. For this reason we 
are not able to compare the results of our investigation to ex- 
tensive theoretical predictions for individual elements in the at- 
mospheres of Ap stars. The standard diff'usion calculations con- 
sider chemical transport processes in the presence of axisym- 
metric magnetic field and thus anticipate abundance inhomo- 
geneities to be axially symmetric themselves (Michaud et al. 
I1981I I. These predictions of the simple diffusion theory are in- 
consistent with the MDI results obtained for HD 24712. We find 
a nearly dipolar magnetic field geometry coexisting with diverse 
abundance distributions, frequently shifted relative to the mag- 
netic field axis. A similar inconsistency between a simple mag- 
netic field topology and complex abundance distributions was 
found in other recent Dl studies of Ap stars (Kochukhov et al. 
120021 |2004bl l. It is possible that diffusion processes are sensitive 
to the small scale magnetic structures that could be resolved only 
with the full Stokes vector MDI, such as the study carried out by 
Kochukhov et al. (2004a) for 53 Cam. On the other hand one 
can also suspect that the current simple diffusion theory lacks 
certain physical processes other than the magnetic field, that can 



influence or even drive chemical spot formation. Recent discov- 
ery of the evolving chemical clouds in the atmosphere of non- 
magnetic chemically peculiar stars (Kochukhov et al. | 2007| l has 
demonstrated that the surface structure formation in the upper 
main sequence stars is not necessarily driven by the magnetic 
field alone. 

Anticipating further development of theoretical aspects, we 
provide new information for the diffusion properties within the 
atmosphere of a cool Ap star, obtain new insights in the atmo- 
spheric structure and the field geometry, the origin and interplay 
of abundance variations, pulsation and magnetic fields in these 
unique stellar laboratories. 
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Table 3. Table of elements and spectral lines used for mapping fromleft to right: element identification, centi-al wavelength, the 
loggf values of these lines, abundance interval and the range in dex, and the solar value as a comparison presented in columns from 
left to right. 



Element 


1 ^ J \ 
A (A) 


loggf 


min, max 


abundance interval 


O 








log(7V/Mot) 


log(Af/Aftot) 


logCW/Mot) 


Mgi 


5528.405 


-0.620 


-6.0, -4.5 


1.5 


-4.51 


Ca n 


5021.138 


-1.456 


-5.2,-5.0 


0.2 


-5.73 


Sc II 


5526.790 


0.024 


-10.1,-8.9 


1.2 


-8.99 


Ti n 


4805.085 


-0.960 


-7.8,-7.1 


0.7 


-7.14 




4911.195 


-0.610 








Cr II 


5510.702 


-2.614 


-6.1,-5.9 


0.2 


-6.40 


Pel 


5005.711 


-0.180 


-5.3,-4.7 


0.6 


-4.59 




5006.117 


-0.638 










5242.490 


-0.967 










5263.300 


-0.879 










5263.394 


-2.774 










5367.470 


0.443 










5383.370 


0.645 










5397.130 


-1.993 










5397.193 


-0.342 








Coi 


4813.449 


-2.121 


-6.4, -5.8 


0.6 


-7.12 




4813.467 


0.050 










4897.179 


-0.764 










5146.161 


-3.660 










5146.749 


-0.320 










5146.823 


-2.393 










5146.998 


-1.286 










5250.000 


0.320 










5378.247 


-0.425 








Nil 


5035.357 


0.290 


-7.6,-6.3 


1.3 


-5.81 




5146.480 


-0.060 








Y n 


5119.112 


-1.360 


-9.0,-8.5 


0.5 


-9.83 




5402.774 


-0.510 








La n 


4804.039 


-1.494 


-9.8,-8.8 


1.0 


-10.91 




4808.996 


-1.494 










4921.776 


-0.699 










5377.052 


-0.316 








Ce n 


5117.169 


-0.050 


-9.3, -8.9 


0.4 


-10.46 




5273.899 


-1.210 










5274.229 


0.150 








Pr II 


5135.140 


-0.131 


-10.0,-9.6 


0.4 


-11.33 


Pr III 


4910.823 


-1.954 


-8.4, -7.7 


0.7 






4929.115 


-2.068 








Nd n 


4811.342 


-1.015 


-9.3,-8.5 


0.8 


-10.59 




5033.507 


-0.470 










5132.328 


-0.710 










5143.337 


-1.570 










5276.869 


-0.393 










5361.467 


-0.482 










5385.888 


-0.860 










5399.099 


-1.419 








Nd m 


4796.499 


-1.657 


-8.0,-6.9 


1.1 


-10.59 




4821.980 


-2.405 










4911.653 


-1.826 










4912.943 


-1.789 










4927.488 


-0.850 








Gd n 


5500.419 


-1.333 


-8.8,-8.5 


0.3 


-10.92 


Tb m 


5505.409 


-0.815 


-9.1,-7.7 


1.4 


-11.76 


Dy n 


4923.167 


-2.384 


-8.8,-8.3 


0.5 


-10.90 
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Fig. 9. Effect of a small rotation of the surface structure of Y, La, and Ce on the Stokes / profiles of the lines of these elements used 
for mapping. Panels show the difference between the original synthetic profiles and spectra computed with the surface spots shifted 
by 0.05 (blue) and 0.10 (red) of the rotation period (18 and 36° on the stellar surface, respectively). The grey bar indicates the 
uncertainty (±cr) of the spectra corresponding to the S /N of individual rotation phases. The difference spectra are offset vertically 
for display purposes. The horizontal scale is arbitrary. 



